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We have used three large N-body simulations to investigate the dependence of dark matter 
halo concentrations on halo mass and redshift in the WMAP year 3 cosmology. The median 
concentration - mass relation is well fit by a power law over the mass range 10'! —10'° A~'Mo 
at redshifts z < 2, for both NFW and Einasto density profiles. We compare our findings to 
recent work based on the Millennium Simulation, which uses a higher value of og, consistent 
with the WMAP year 1 cosmology. We find that for z = 0, halo concentrations are reduced 
by up to 34 per cent at 10''h~+Mo, and 19 per cent at 10'4h71Mo. 


1 Introduction 


It has been argued! that the concentration of a halo reflects the background density of the 
Universe at its formation time, and as small objects form first in a hierarchical universe, lower 
mass halos will be more concentrated than massive ones. 

We have used three DM only simulations, run with GADGET2?, Each simulation contains 
512? DM particles but with a progressively larger comoving box size: 25, 100, and 400h~!Mpc 
for runs L025, L100 and L400, respectively. We have used the 3 year Wilkinson Microwave 
Anisotropy Probe (WMAP) cosmology®, hereafter WMAP3, and compared with the Millennium 
Simulation* (MS) c(M) relation® for the first year WMAP release’ henceforth WMAP1; where 
og is about 23 per cent larger, therefore one can expect a given mass halo to form at earlier 
times and hence with higher concentration. Please refer to Duffy et al? for further details. 


1.1 Density Profiles 


The Navarro-Frenk-White | (NFW) density profile is given by 
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where 6, is a characteristic density contrast and rs is a scale radius. The concentration parameter 
is defined as c200 = r200/rs. NFW profiles were fit using the two parameters rs and ôe. Other 
concentrations, e.g. Cyir and Cmean, can also be defined for different overdensity values. 


The Einasto profile is a rolling power-law first introduced to describe the distribution of old 
stars in the Milky Way ®. It takes the form 


2 
In(p/p_2) = =z [(r/r-2)* — 1] , (2) 
where p_2 is the density at r- the radius where ge = —2 and a is the rate at which the 


profile turns over. An analagous concentration to the NFW case can be defined using rs. 
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Figure 1: Concentration-mass relations for z = 0 (top), 1 (middle) and 2 (bottom) using NFW density profiles. 

Data points correspond to median values and error bars to quartiles. Only bins containing at least 5 haloes are 

shown. The solid lines show the best-fit power-law relation. The errors on the best-fit parameters, given in the 

legend, are determined by bootstrap sampling the haloes and correspond to 68.2 per cent confidence limits. The 
dashed line in the left panel shows the best-fit power-law relation from Neto et al. 


Table 1: Best-fit parameters for the median c(M,z) relation using Mpivot = 2 X 1012 h7'Mo for three halo 
definitions using two different density profiles. The errors correspond to 1o confidence intervals and have been 
determined by bootstrap resampling the halos. 


Einasto 


2 Results 


The data, which spans four orders of magnitude in mass and a redshift range z = 0 — 2 are well 
described by a function of the form 


c = A(M/Mpivot)?(1 + 2)”. (3) 


The parameter values and errors are given in Table 1. 
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